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ALMOST-SURE GROWTH RATE OF GENERALIZED RANDOM FIBONACCI 

SEQUENCES 

ELISE JANVRESSE, BENOIT RITTAUD, THIERRY DE LA RUE 



Abstract. We study the generalized random Fibonacci sequences defined by their first non- 
negative terms and for n > 1, -Fn+2 = A_F„+i ± Fn (linear case) and Fn+2 = \^Fn+\ ± Fn\ 
(non-linear case), where each ± sign is independent and either + with probability p or — with 
jH ' probability 1 — p (0 < p < 1). Our main result is that, when A is of the form Aj; = 2cos(7r/fc) 

for some integer fc > 3, the exponential growth of Fn for < p < 1, and of Fn for 1/fc < p < 1, 
is almost surely positive and given by 



/■oo 

Jo 



where p is an explicit function of p depending on the case we consider, taking values in [0, 1], 
and ^j. p is an explicit probability distribution on R-|_ defined inductively on generalized Stern- 

P^ ^ Brocot intervals. We also provide an integral formula for < p < 1 in the easier case A > 2. 

P^ . Finally, we study the variations of the exponent as a function of p. 

1. Introduction 

Random Fibonacci sequences have been defined by Viswanath by Fi = i^2 = 1 and the random 
recurrence Fn+2 = -Fri+i ± Fm where the ± sign is given by tossing a balanced coin. In [11], he 
^T . proved that 

^; a/I^ — ^1.13198824... a.s. 

CO ' and the logarithm of the limit is given by an integral expression involving a measure defined on 

Cn \ Stern-Brocot intervals. Rittaud ^ studied the exponential growth of E(|F„|): it is given by an 

^^ ■ explicit algebraic number of degree 3, which turns out to be strictly larger than the almost-sure 

^D \ exponential growth obtained by Viswanath. In [S] , Viswanath's result has been generalized to the 

QP ■ case of an unbalanced coin and to the so-called non-linear case Fn+2 = \Fn+i ± ^n|- Observe that 

this latter case reduces to the linear recurrence when the ± sign is given by tossing a balanced 
coin. 

A further generalization consists in fixing two real numbers, A and /3, and considering the 
5^ ■ recurrence relation Fn+2 = AF„+i ± (iFn (or Fn+2 = |Ai^n+i ± /3-Fn|), where the ± sign is chosen 

5^ \ by tossing a balanced (or unbalanced) coin. By considering the modified sequence G„ := i^„//3"/^, 

which satisfies Gn+2 — -y^Gn+i ± G„, we can always reduce to the case /3 = 1. The purpose of 
this article is thus to generalize the results presented in [5] on the almost-sure exponential growth 
to random Fibonacci sequences with a multiplicative coefficient: (i^n)n>i and {Fn)n>i, defined 
inductively by their first two positive terms Fi ~ Fi = a, F2 = F2 = b and for all n > 1, 

(1) Fn+2 = Ai^„+i ± Fn (linear case), 

(2) Fn+2 = |AF„+i ± Fn\ (non- linear case), 

where each ± sign is independent and either + with probability p or — with probability 1 — p 
(0 < p < 1). We are not yet able to solve this problem in full generality. If A > 2, the linear 
and non-linear cases are essentially the same, and the study of the almost-sure growth rate can 
easily be handled (Theorem II. 3p . The situation A < 2 is much more difficult. However, the 
method developed in [Hj can be extended in a surprisingly elegant way to a countable family of 
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A's, namely when A is of the form A^, = 2cos(7r/fc) for some integer fc > 3. The simplest case 
A3 = 1 corresponds to classical random Fibonacci sequences studied in [5]. The link made in [S] 
and [8] between random Fibonacci sequences and continued fraction expansion remains valid for 
Afc = 2 cos(7r/fc) and corresponds to so-called Rosen continued fractions, a notion introduced by 
Rosen in [9] . These values Afc are the only ones strictly smaller than 2 for which the group (called 
Hecke group) of transformations of the hyperbolic half plane H^ generated by the transformations 
z I — > ^1/2 and z 1 — > z + A is discrete. 

In the linear case, the random Fibonacci sequence is given by a product of random i.i.d. ma- 
trices, and the classical way to investigate the exponential growth is to apply Furstenberg's for- 
mula [3j. This is the method used by Viswanath, and the difficulty lies in the determination of 
Furstenberg's invariant measure. In the non- linear case, the involved matrices are no more i.i.d., 
and the standard theory does not apply. Our argument is completely different and relies on some 
reduction process which will be developed in details in the linear case. Surprisingly, our method 
works easier in the non-linear case, for which we only outline the main steps. 

Our main results are the following. 

Theorem 1.1. Let A = A^ = 2cos(7r/fc), for some integer fc > 3. 

For any p € [0,1], there exists an explicit probability distribution vj^ p on R^_ defined inductively 
on generalized Stern-Brocot intervals (see Section 1X21 and Figurc[Tjl, which gives the exponential 
growth of random Fibonacci sequences: 

• Linear case: Fix Fi > and F2 > 0. For p = 0, the sequence {\Fn\) is periodic with 
period fc. For any p g]0, 1], 

1 f°° 

-log|F„| > jp_Xk = / logxdi^k.p{x) > 







almost-surely, where 



and PR is the unique positive solution of 

fe-i 
1 P- ^ ] ^1 



• 



p+ {1 — p)x ^ 

Non-linear case: For p G]l/fc, 1] and any choice of Fi > and F2 > 0, 

1 ~ ^ f°° 

-logFn > 7p.A^ = / log xdi^k,p{x) > 

n n^oo ' Jg 

almost-surely, where 

p:= "-^I-PR 
and pr is, for p < 1, the unique positive solution of 



k-l 
1-—^ ] =1-,. 

(1 -p) +px^ 



(Forp= l,pfl = 1.; 



The behavior of (Fn) when p < 1/fc strongly depends on the choice of the initial values. This 
phenomenon was not perceived in [S], in which the initial values were set to i^i = F2 = 1. However, 
we have the general result: 

Theorem 1.2. Let A = Afc = 2cos(7r/fc), for some integer fc > 3. In the non-linear case, for 
<p < 1/fc, there exists almost-surely a bounded subsequence {Fn ) of{Fn) with density [l — kp). 



The bounded subsequence in Theorem 11.21 satisfies Fn^^-^ = |AF„^ — Fn^_-^\ for any j, which 
corresponds to the non- linear case for p = 0. We therefore concentrate on this case in Section W?I 
and provide necessary and sufficient conditions for {Fn) to be ultimately periodic (see Proposi- 
tion lG.Sp . Moreover, we prove that Fn may decrease exponentially fast to 0, but that the exponent 
depends on the ratio Fq/Fi. 
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Figure 1. The measure v^^p on generalized Stern-Brocot intervals of rank 1 and 
2 in the case fc = 4 (Afe = v^)- The normalizing constant Z is given by \ + p + p^. 
The endpoints of the intervals are specified by their \/2-continued fraction expan- 



sion. 



The critical value 1/fc in the non-linear case is to be compared with the results obtained in 
the study of E[i^„] (see [4]): it is proved that E[i^„] increases exponentially fast as soon as 
p > (2 - Afe)/4. 

When A > 2, the linear case and the non-linear case are essentially the same. The study of the 
exponential growth of the sequence (fn) is much simpler, and we obtain the following result. 

Theorem 1.3. Let A > 2 and < p < 1. For any choice of Fi > and F2 > 0, 
1 



■log|i^„| 



7p,A 



log xdfip^\{x) > 



where jip^x is an exphcit probabihty measure supported on \B, A -I- -i] , with B 
(see Section 13 and Figure\^. 



A + VA2 - 4 



Road map. The detailed proof of Theorem 11.11 in the linear case is given in Sections [^HH Sec- 
tion [2] explains the reduction process on which our method relies. In Section [31 we introduce the 
generalized Stern-Brocot intervals in connection with the expansion of real numbers in Rosen con- 
tinued fractions, which enables us to study the reduced sequence associated to (Fn). In ScctionlH 
we come back to the original sequence (Fn), and, using a coupling argument, we prove that its 
exponential growth is given by the integral formula. Then we prove the positivity of the integral 
in Section [5] 

The proof for the non-linear case, p > 1/fc, works with the same arguments (in fact it is even 
easier) , and the minor changes are given at the beginning of Section [6] The end of this section is 
devoted to the proof of Theorem 11.21 

The proof of Theorem 1 1.31 (for A > 2) is given in Section [71 

In Section FS.ll we study the variations of jp^\ and jp,\ with p. Conjectures concerning variations 
with A are given in Section \K% 

Connections with Embree-Trefethen's paper [2], who study a slight modification of our linear 
random Fibonacci sequences when p = 1/2, arc discussed in Section [9l 



2. Reduction: The linear case 

The sequence (-Fn)ri>i can be coded by a sequence (X„)„>3 of i.i.d. random variables taking 
values in the alphabet {R,L} with probability {p, 1 — p). Each R corresponds to choosing the + 
sign and each L corresponds to choosing the — sign, so that both can be interpreted as the right 
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multiplication of (i^„_i, F^) by one of the following niatriees: 

(3, ..(; -) .„a «.(; I 

Aecording to the context, we will interpret any finite sequence of i?'s and i's as the corresponding 
product of matrices. Therefore, for all n > 3, 

{Fn-l, Fn) ~ (Fi, ^2)^3 . . . Xn- 

Our method relies on a reduction process of the sequence (X„) based on some relations satisfied 
by the matrices R and L. Recalling the definition of A = 2cos(7r/fc), we can write the matrix L 
as the product P^^DP, where 

D := \ „ _,,./t , P := L _,^ii.\, and P '- = 



e-"/V ' • Vl e-^V ' 2isin(7r/A;) V 1 "1 

As a consequence, we get that for any integer j, 



i / — sm " , ' 



sm(7r/K) \ sm ^ sm 




and 



(5) RV 



1 / • 77r 

i / sm ''t sm 



JTT ■ (j + 1) 



k k 



sin(7r/fc) Isinii!^ sin ^^ 



In particular, for j = fc — 1 we get the following relations satisfied by R and L, on which is based 
our reduction process: 

(6) RL''-^ = \\ ^V RL^-^R = -L and RL^-^L=-R. 

Moreover, L^ = — Id. 

We deduce from ^ that, in products of i?'s and L's, we can suppress all patterns RL*'~^ 
provided we flip the next letter. This will only affect the sign of the resulting matrix. 

To formalize the reduction process, we associate to each finite sequence x — x^ . . .Xn G 
{i?, L}"~^ a (generally) shorter word Red(a:) = y3---yj by the following induction. If n = 3, 
2/3 = 2:3. If n > 3, Red(a;3 . . . x„) is deduced from Red(a:3 . . . x„_i) in two steps. 

Step 1: Add one letter (i? or L, see below) to the end of Red(x3 . . . Xn-i)- 

Step 2: If the new word ends with the sufRx RL''~^, remove this suffix. 

The letter which is added in step 1 depends on what happened when constructing Red (0:3 . . . x„_i): 

• If Red(a;3 . . . x„_i) was simply obtained by appending one letter, we add Xn to the end of 
Red(a;3 . . .a;„_i). 

• Otherwise, we had removed the suffix RL^^^ when constructing Red(a;3 . . . a;„„i); we then 
add T^ to the end of Red(a;3 . . . x„_i), where R:= L and L := R. 

Example: Let x = RLRLLLRLL and fc = 4. Then, the reduced sequence is given by Red(x) = R. 

Observe that by construction, Red(x) never contains the pattern RL''~^ . Let us introduce the 
reduced random Fibonacci sequence {F^) defined by 

(K-i,J^;) := (Pi,P2)Red(X3...X„). 

Note that we have Fn = ^F^ for all n. From now on, we will therefore concentrate our study 
on the reduced sequence Red(X3 . . . Ar„). We will denote its length by j{n) and its last letter by 
Y{n). 

The proof of Lemma 2.1 in |^ can be directly adapted to prove the following lemma. 

Lemma 2.1. We denote by \W\r the number of R's in the word W. We have 

(7) \Rcd{X3...Xn)\R >+^ a.s. 
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In particular, the length j{n) of Rcd{X^ . . . X„) satisfies 

j{n) > +00. a.s. 

n — ^00 

2.1. Survival probability of an R. We say that the last letter of Red(X3 . . .X„) survives if, 
for all m > n, j{rn) > j{n). In other words, this letter survives if it is never removed during 
the subsequent steps of the reduction. By construction, the survival of the last letter Y{n) of 
Red(X3 . . . Xn) only depends on its own value and the future X„+i, X,i+2 ■ • ■• Let 



Pb, := P( y{n) survives y{n) = R has been appended at time n\ . 



A consequence of Lemma 12711 is that pR > 0. We now want to express pR as a function of p. 

Observe that y{n) = R survives if and only if, after the subsequent steps of the reduction, it 
is followed by U R where < j < fc — 2, and the latter R survives. Recall that the probability of 
appending an R after a deletion of the pattern RL^"^ is 1 — p, whereas it is equal to p if it does 
not follow a deletion. Assume that Y{n) = R has been appended at time n. We want to compute 
the probability for this R to survive and to be followed by L^ R {0 < j < k — 2) after the reduction. 
This happens with probability 



p, := F 



/ -" \ 

/ deletions survives \ 

i? be followed by ([R^l] ... ([R^L] [R^ R 



\ j times / 



J 



(1 -p) +p^(i -pj,yii-pY~'p pY,i^-PRYii~pYpR 



1- 



e>i I i>Q 

PPR \ ' PPR 



P + (1 - p)pR J P + (1 - P)PR ' 

Writing pR ~ X) i=o Pj^ ^® S,^^ that pR is a solution of the equation 

(8) g(x) = 0, where g{x) := 1 ^^-^ - (1 - x)^/'^''-^^ . 

Observe that g{0) = 0, and that g is strictly convex. Therefore there exists at most one a; > 
satisfying g{x) — 0, and it follows that pR is the unique positive solution of ([8]). 



2.2. Distribution law of surviving letters. A consequence of Lemma [2T] is that the sequence 
of surviving letters 

(5,)j>3= lim Rcd(X3...X„) 

n — 'oo 

is well defined and can be written as the concatenation of a certain number s > of starting L's, 
followed by infinitely many blocks: 

Si 02 . . . ^ L B1B2 . . . 

where s > and, for all i > 1, Bg E {R, RL, . . . , RL''^^}. This block decomposition will play a 
central role in our analysis. 

We deduce from Section ETT] the probability distribution of this sequence of blocks: 

Lemma 2.2. Tiic blocks {Bi)i>i are i.i.d. with common distribution law Pp defined as follows 
(9) Pp(Si = RL^) := J^l , < J < fc - 2, 

where p '.= 1 -, -^ — and pR is the unique positive solution of dSl). 

p + (1 - pjpR 
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In [S], where the case fc = 3 was studied, we used the parameter a ~ 1/(1 + p) instead of p. 

Observe that p = ( (1 — p) + P^gyi{l — PnYi^ ^ pY^^p) can be interpreted as the probabihty 
that the sequence of surviving letters starts with an L. Since an R does not survive if it is followed 
by fc — 1 L's, this explains why the probability 1 — p/j that an R does not survive is equal to p^~^ . 

Proof. Observe that the event En :="y(n) = R has been appended at time n and survives" is the 
intersection of the two events "F (n) = R has been appended at time n" , which is measurable with 
respect to cr(Xj, i < n), and "If Y{n) = R has been appended at time n, then this R survives", 
which is measurable with respect to cr(Xi, i > n). It follows that, conditioned on E^ o'(Xj, i < n) 
and a{Xi, i > n) remain independent. Thus the blocks in the sequence of surviving letters appear 
independently, and their distribution is given by 



Pp(Bi = RV) = ^ = —^ , {)<]<k-2. 



3. Rosen continued fractions and generalized Stern-Brocot intervals 



n 



3.1. The quotient Markov chain. For £ > 1, let us denote by ng the time when the £-th 
surviving R is appended, and set 

n«+l-2 

Qi is the quotient of the last two terms once the £-th definitive block of the reduced sequence has 
been written. Observe that the right-product action of blocks B € {i?, RL, . . . , RL^~^} acts on 
the quotient -F^/F^.i in the following way: For < j < fc — 2, for any (a, 5) G R* x R, if we set 
(a', 6') := {a,b)RU, then 

a' \a 

where fo{q) := A + 1/q and f{q) := A — 1/q. For short, we will denote by fj the function /-' o /o. 
Observe that fj is an homographic function associated to the matrix RL^ in the following way: 

To the matrix „ J corresponds the homographic function q i-^ ^J^ "^ . 

It follows from Lemma 12.21 that {Qe)e>i is a real- valued Markov chain with probability transi- 
tions 

P (Q,+i = fj{q)\Qi ^q) = — fc^ , < .7 < fc - 2. 

Z^rn=0 P 

3.2. Generalized Stern-Brocot intervals and the measure Vk.p- Let us define subintervals 
of R: for < j' < fc — 2, set Ij := /j([0, 4-oo]). These intervals are of the form 

I'j = [bj+i^bj], where 6o = +oo, &i = A = /o(+oo) = /i(0), bj+i = f{hj) = fji+oo) = fj+i{0). 

Observe that bk~i = ./fe-i(O) = since RL''^^ ~ [n i )• Therefore, iIj)o<j<k-2 is a subdivi- 
sion of [0, -l-oo]. 

More generally, we set 

Iji.h.-.n ■= fji ° /j2 ° ■ ■ • ° /j>([0, +oo]), V(ji, j2, . . . , jf) e {0, . . . , fc - 2Y. 
For any £ > I, this gives a subdivision J^{i) of [0, +oo] since 

fc-2 
3e=0 

When fc = 3 (A = 1), this procedure provides subdivisions of [0, +oo] into Stern-Brocot intervals. 
Lemma 3.1. The a-algebra generated by J^{£) increases to the Borel a-algebra on R+. 
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Wc postpone the proof of this lemma to the next section. 

Observe that for any q G R+. W{Qi G Iji.j2 jilQo = q) = f^k-2 „,w - Therefore, the proba- 

bihty measure i>k,p on R-|- defined by 

„ilH hji 

^k,p\J-ji,J2,...,ji) • = 



(Eo p")' 

is invariant for the Markov chain [Qi). The fact that Vk,p is the unique invariant probabihty for 
this Markov chain comes from the foUowing lemma. 

Lemma 3.2. There exists almost surely L+ > such that for all £ > L+, Q( > 0. 

Proof. For any q G R \ {0}, cither fo{q) > 0, or fi{q) = A - l/faiq) > 0. Hence, for any i>0, 

F{Qe+i > 0\Qi ^q)> — ^^ . 

It follows that F(V£ > 0, Qe < 0) ~ 0, and since Q^ > =^ Qe+i > 0, the lemma is proved. D 

To a given finite sequence of blocks {RL^'^), . . . , {RU^), we associate the generalized Stern- 
Brocot interval Iji,j2,---,3i- If ^^ extend the sequence of blocks leftwards, we get smaller and 
smaller intervals. Adding infinitely many blocks, we get in the limit a single point corresponding 
to the intersection of the intervals, which follows the law i^t.p. 

3.3. Link with Rosen continued fractions. Recall (see [9j) that, since 1 < A < 2, any real 

number q can be written as 

1 
q — aoX H 

fliA H 



1 



a„A+ 



where (a„)„>o is a finite or infinite sequence, with a„ G Z \ {0} for n > 1. This expression will 
be denoted by [ao, . . . , a„, . . .]a. It is called a X-Rosen continued fraction expansion of q, and is 
not unique in general. When A = 1 (i.e. for fc = 3), wc recover generalized continued fraction 
expansion in which partial quotients are positive or negative integers. 

Observe that the function fj are easily expressed in terms of Rosen continued fraction expansion. 
The Rosen continued fraction expansion of fj (g) is the concatenation of {j + 1) alternated ±1 with 
the expansion of ±g according to the parity of j: 



(10) fj{[ao,...,an,...U) = < 



[1,-1,1,.. .,l,ao,.. .,a„,.. .]a if j is even 

ij+i) terms 
[1,-1,1,.. .,-1,-ao,. .. ,-a„,. ..]a if j is odd. 



(j+i) terms 

For any ^ > 1, let £"{£) be the set of endpoints of the subdivision y{i)- The finite elements of 
(f(l) can be written as 

j terms 

In particular for j = A: — 1 wc get a finite expansion of bk-i = 0. Moreover, by pop . 

bo = /o(0) = c» = [1, 1,-1,1 , ...,±1 ]a. 

fc-1 terms 

Iterating (fTO|) . we see that for all £ > 1, the elements of (?(£) can be written as a finite A- Rosen 
continued fraction with coefficients in { — 1, 1}. 

Proposition 3.3. The set lj£>i '^e of all endpoints of generalized Stern-Brocot intervals is the set 
of all nonnegative real numbers admitting a finite X-Rosen continued fraction expansion. 
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The proof uses the two following lemmas 
Lemma 3.4. 



Proof. From ([5]), we get 






-k-2 /A A , . _ 1 



^^ =1 n ' hence fk-2ix) 



1 0/' ^---v-^ ;^^^- 



Therefore. fk-2{^/x) = 1/ fo{x) and the statement is true for j = 0. Assume now that the result 
is true for j > 0. We have 

so the result is proved by induction. D 

Lemma 3.5. For any i > 1, the set S'{£) of endpoints of the subdivision J^{£) is invariant by 
X i-^ 1/x. Moreover, the largest finite element of £"{£) is iX and the smallest positive one is 1/£X. 

Proof. Recall that the elements of 1^(1) are of the form bj = /,_i(oo) = /j(0), and the largest 
finite endpoint is bi = A. Hence, the result for ^ = 1 is a direct consequence of Lemma 13.41 

Assume now that the result is true for i > 1. Consider b G S'{i + 1) \ 'S'{i). There exists 
< J < k — 2 and 6' G £"{£) such that b = fj{b'). Since 1/6' is also in (^{i), we see from Lemma [5^ 
that 1/b = fk-2-ji^/b') S £"{£ + 1). Hence S'{£ + 1) is invariant by a; f-> 1/x. Now, since /o is 
decreasing, the largest finite endpoint of £"{£ + 1) is /o(l/^A) = {£+ 1)A, and the smallest positive 
endpoint of y{£ + 1) is l/{£ + 1)A. D 



Proof of Proposition \S.3[ The set of nonnegative real numbers admitting a finite A- Rosen contin- 
ued fraction expansion is the smallest subset of IR.)- containing which is invariant under x i— > 1/x 
and X i—f X + X. By Lemma 13.51 the set 1J^>]^ (o/ is invariant under x i— > 1/x. Moreover, it is also 
invariant by a; i— + fk-2{x) = l/(a; + A), and contains 6fe_i =0. D 

Remark 3.6. The preceding proposition generalizes the well-known fact that the endpoints of 
Stern-Brocot intervals are the rational numbers, that is real numbers admitting a finite continued 
fraction expansion. 

Proof of Lemma \3.1[ This is a direct consequence of Proposition 13.31 and the fact that the set of 
numbers admitting a finite A-Rosen continued fraction expansion is dense in R for any A < 2 (see 

m- □ 

4. Coupling with a two-sided stationary process 

If |F„+i/F„| was a stationary sequence with distribution h'k,p, then a direct application of the 
ergodic theorem would give the convergence stated in Theorem 11.11 The purpose of this section 
is to prove via a coupling argument that everything goes as if it was the case. For this, we embed 
the sequence (X„)„>3 in a doubly- infinite i.i.d. sequence (X*)„ga with X„ = X* for all n > 3. 
We define the reduction of {X*)^co<j<n, which gives a left-infinite sequence of i.i.d. blocks, and 
denote by g* the corresponding limit point, which follows the law Vk.p- We will see that for n large 
enough, the last £ blocks of Red(Ar3 . . . A"„) and Red((A'*)_oo<j<n) are the same. Therefore, the 
quotient (/„ := F^/F^_i is well- approximated by g* , and an application of the ergodic theorem to 
g* will give the announced result. 

4.1. Reduction of a left-infinite sequence. We will define the reduction of a left-infinite i.i.d. 
sequence (X*)^^ by considering the successive reduced sequence Red(Ar^„ . . .Xq). 

Proposition 4.1. For all £ > 1, there exists almost surely N{£) such that the last £ blocks of 
Red(X:i„ . . . Xq) are the same for any n > N{£). 
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This allows us to define almost surely the reduction of a left-infinite i.i.d. sequence (X*)'^^^ as 
the left-infinite sequence of blocks obtained in the limit of Red(Xl„ . . . Xq) as n — > oo. 

Let us call excursion any finite sequence wi . . . Wm of R's and L's such that Red(it;i . . . Wm) = 0- 
We say that a sequence is proper if its reduction process does not end with a deletion. This means 
that the next letter is not flipped during the reduction. 

The proof of the proposition will be derived from the following lemmas. 

Lemma 4.2. If there exists n > such that Xl„ . . . X';^-^ is not proper, then Xq is preceded by a 
unique excursion. 

Proof. Wc first prove that an excursion can never be a suffix of a strictly larger excursion. Let 
W = WiRW be an excursion, with RW another excursion. Then WL ~ WiRW L ~ ±R and 

RW'L = ±i?, which implies that Wi = ±Id. It follows that Red(W^i) = [} i*-, ) • Observe 

that Red(iyi) cannot start with L's since Ked{WiRW') = 0. Therefore, it is a concatenation of 
s blocks, corresponding to some function fj-^ o ■ ■ ■ fj^ which cannot be a; i-^ ±x unless s = 0. But 
s = means that Red(Tyi) = 0, so Red(VF) == Rcd{LW') = 0, which is impossible. 

Observe first that, if Xq is not flipped during the reduction of X'^,_^. . . .Xq but is flipped 
during the reduction of X*„...Xq, then Xl„ is an R which is removed during the reduction 
process of X^„ . . .Xq. In particular, this is true if we choose n to be the smallest integer such 
that Xq is flipped during the reduction of Xl„ . . -^o- Therefore there exists < j < n such 
that XZn ■ ■ ■ ^1( .Li-i is an excursion. If j = we are done; otherwise the same observation proves 
that X'^, is an L which is flipped during the reduction process of X^„ . . . XI -. Therefore, Xq is 
flipped during the reduction oi RX'^,-^^ . . . Xq, but not during the reduction oi X'^^ . . . Xq for any 
i < j— 1. Iterating the same argument finitely many times proves that Red(Xl„ . . . X^i) = 0. D 



Lemma 4.3. 



"S^ ¥{w) < 1. 



w excursions 



Proof. Xq is an R which does not survive during the reduction process if and only if it is the 
beginning of an excursion. By considering the longest such excursion, we get 

p{^-Pr)= X! ^(w) {l-p)PR+P 
Hence, 



w excursions 



(11) E ^H = 71 



< 1. 



w excursions 



P(l - Pb) 

{1-P)PR+P 

n 



We deduce from the two preceding lemmas: 



Corollary 4.4. There is a positive probabihty that for all n > the sequence X^„ . . .Xt^ be 
proper. 



Proof of Proposition \4-l\ We deduce from Corollary 14.41 that with probability 1 there exist infin- 
itely many j's such that 

• XZ j is an R which survives in the reduction of X^ ■ . . . Xq ; 

• X^,j . . . -'^Ij-i is proper for all n > j. 

For such j, the contribution of XI ■ . . . Xq to Red(X^„ . . . Xq) is the same for any n > j. D 

The same argument allows us to define almost surely Red((X*)"g^) for all n E TL, which is 
a left-infinite sequence of blocks. Observe that we can associate to each letter of this sequence 
of blocks the time t < n at which it was appended. We number the blocks by defining Bq as 
the rightmost block whose initial R was appended at some time < < 0. For n > 0, we have 
Rcd((X*)!!^) = . . . BI^B^B'!^ . . . B2(„) where < L{n) < n. The random number L{n) evolves 
in the same way as the number of R's in Red(X3 . . . X„). By Lemma |2.1[ L(n) —f -l-oo as n — > (X) 
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almost surely. As a consequence, for any j G Z the block B" is well-defined and constant for 
all large enough n. We denote by Bj the limit of B". The concatenation of these blocks can be 
viewed as the reduction of the whole sequence (X*)^^. The same arguments as those given in 
Section [2] prove that the blocks Bj are i.i.d. with common distribution law P^. 

It is remarkable that the same result holds if we consider only the blocks in the reduction of 

(^*)°-oo- 

Proposition 4.5. The sequence Red{{X*)^^) is a left-inEnite concatenation of i.i.d. blocks with 
common distribution law P^. 

Proof. Observe that Red((X*)'L„^) = Red{{X*)'l^) where L < is the (random) index of the 
last letter not removed in the reduction process of (X*)'^^. For any £ < 0, we have i = ^ if and 
only if {X*Y_^ is proper and {X*)'^_^_^ is an excursion. For any bounded measurable function /. 
since ]E,[f {Rcd{{X*Y_^)) | {X*Y_^ is proper] does not depend on £, we have 

E[/(Rcd((X*)°^)] 

= ^P(i = £)E[/(Red((X*)i^))|i = ^] 

i 

= ^F{L = e) E[/(Rcd((X*)i^)) I (X*)i^ is proper, (X*)^_^i is an excursion] 

e 

= Y.^{L = i) E[/(Red((X*)i^)) I iX*Y-^ is proper] 

= E[/(Red((X*)«_^)) I {X*Y.^ is proper]. 

This also implies that the law of Red((X*)^o^) is neither changed when conditioned on the fact 
that {X*Y_^ is not proper. 

Assume that (X*)°^ is proper. The fact that the blocks of Red((X*)^^) will not be sub- 
sequently modified in the reduction process of (X*)^^ only depends on (X*)^. Therefore, 
E[/(Red((A:*)0 ^)) I iX*Y_^ is proper] is equal to 

E[/(Red((X*)°^)) I iX*Y_^ is proper and blocks of Red((X*)°^) are definitive]. 

The same equality holds if we replace "proper" with "not proper". Hence, the law of Red((A'*)" q^) 
is the same as the law of Red((X*)'io^) conditioned on the fact that blocks of Rcd((X*)° g^) are 
definitive. But we know that definitive blocks are i.i.d. with common distribution law Pp. D 

4.2. Quotient associated to a left-infinite sequence. Let n be a fixed integer. For to > 0, we 
decompose Red{{X*)n-m<i<n) into blocks B^, . . . ,Bi = [RU'^), . . . , {RU^), to which we associate 
the generalized Stern-Brocot interval Iji,j2,---,jf If we let to go to infinity, the preceding section 
shows that this sequence of intervals converges almost surely to a point g* . By Proposition 14.51 
g* follows the law Vk,p- 

Since (g*) is an ergodic stationary process, and log(-) is in L^{i'k.p); the crgodic theorem implies 

1 ^ f 

(12) T^yi^'^Sqn'T: ' ^ogqdvk.piq) almost surely. 

The last step in the proof of the main theorem is to compare the quotient g,i = F[^/Fj[_i with g* . 

Proposition 4.6. 

1 ^ 
— V I log g,* - log I g„ 1 1 > almost surely. 

n— 3 

We call extremal the leftmost and rightmost intervals of J^(£). 

Lemma 4.7. 

Si :— sup sup I logg* — logg| > 

Inot extremal 



ALMOST-SURE GROWTH RATE OF GENERALIZED RANDOM FIBONACCI SEQUENCES 11 

Proof. Fix £ > 0; and choose an integer M > l/e. By Lemma 13. 1[ since log(-) is uniformly 
continuous on [1/il/A, MX], we have for i large enough 

sup sup I logg* — log (/I < e. 

/C[1/MA,MA] 

If / e ^(^) is a non-extremal interval included in [0, l/MX] or in [MA, +00], there exists an integer 
j e [MJ] such that / C [l/(j + 1)A, 1/jA] or / C [{j + 1)A, jA]. Hence, 

sup |logg* -log^l < log (^-^j < log f 1 + -^] <£• 

n 

Proof of Proposition \4-('] For any j G Z, we define the following event Ej: 

• X* is an R which survives in the reduction of (X*)i>j; 

• X* . . . X*_^ is proper for all i < j. 

Observe that if Ej holds for some j > 3, then for all n> j, 

Red(X3...X„) = Red(X3...Xj-i) Red(X,...X„) 

and Red{{X*r_^) ^ Red{{X*yS^) Rcd{X* . . .X:). 

Hence, since Xj . . . Xn = X* . . . X*, they give rise in both reductions to the same blocks, the first 
one being definitive. Since each Ej holds with the same positive probability, the ergodic theorem 
yields 



1 " 
(13) - V Ifi. > FfE's) > almost surely. 



n 

hence the number of definitive blocks of Red(X3 . . . Xn) and of Red((X*)"o^) which coincide grows 
almost surely linearly with n as n goes to cx) (these definitive blocks may be followed by some 
additional blocks which also coincide). 

Recall the definition of L+ given in Lemma \J^ and observe that for n > n^j^, g„ > 0. Observe 
also that, by definition of Iji,j2,...,ji^ if 9 £^nd q* are two positive real numbers, fj-^ o fj^o ■ ■ ■ o fji{q) 
and /jj o /j2 o • • • o fji{q*) belong to the same interval of J'il). 

From p3)) . we deduce that, almost surely, for n large enough, at least L_|_ + ^Jn definitive blocks 
of RcA{Xz ■ . ■ Xn) and of Red((X*)"3^) coincide (possibly followed by some additional blocks 
which also coincide). This ensures that qn and g* belong to the same interval of the subdivision 
J{^). 

By Lemma |4.7[ it remains to check that, almost surely, there exist only finitely many n's such 
that g* belongs to an extremal interval of the subdivision ^(^/n). But this is a direct application 
of Borel-Cantelli Lemma, observing that the measure v^^p of an extremal interval of J" it) decreases 
exponentially fast with i. D 

We now conclude the section by the proof of the convergence to the integral given in Theo- 
rem ll.ll linear case: Since Fn = ^Fn^ we can write n~^ log |F„| as 



1 1 " 1 " 

log 1^2 1 + - E i°g 9.* + - E (i°g 1*1- i°g «;) 



n ■^ — ' ■' n 

and the convergence follows using Proposition 14.61 and p^ . 

5. POSITIVITY OF THE INTEGRAL 

We now turn to the proof of the positivity of jp,\^ . It relies on the following lemma, whose 
proof is postponed. 

Lemma 5.1. Fix < p < 1. For any t > 0, 

(14) At := iyk,p ([i, 00)) - i^k,p ([0, l/t]) > 0. 

Moreover, there exists t > 1 such that the above inequahty is strict. 
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Using Fubini's theorem, we obtain that Jp^\^. is equal to 

/ \ogxdiyk,pix) = / \ogxdvk.pix)- \og{l / x) dvk.p{x) 
Jo Ji Ja 

/•oo 

Vk,p{[e^,oo))du~ I i/fc,p([0,e""])du 

"'0 



which is positive if < p < 1 by Lemma 15.1 1 Thus, ^p.x^, > for any p > 0. This ends the proof 
of Thcorcmll.il hnear case. 



Proof of Lemma I5.il By Lemma 13.11 it is enough to prove the lemma when t is the endpoint of 
an interval of the subdivision J^{i). This is done by induction on £. Obviously, Aq = Aqo = 0. 
When £ = 1 and £ = 2, if i 7^ 0, 00, it can be written as fjibi) for < j < fc — 2 and < i < fc — 2, 
and we get 1/t = fk-2-j{bk-i-i) (see Lemma [3^ . Setting Z :— J2s=o /°'^ ^'^ have 



s^O s— s — s— i 



Therefore, 



i/fc,p([t,oo))-i/fc,p([0,l/t]) 

^zz^z\^z z ^ z 

s— s— i \s — A;— 1— J s— 

Since i < fc - 2, we have p*+J - p'^-^-j > pk-i-j^^^ij _ ^y Moreover, E^^o"* ^ ^ 1- Thus, 

ZA* > ^p^ (1 - p'^-l-J) - p'^-2-,(i _ ^2j)^ 

Observe that (1 - p^-^-') = (1 - p) Y!Zt' P' and that 1 - p^J = (1 + pJ)(l - p) X;lzS p"- Hence, 

^^-2-^"(l + p^)), 



j-l /k-2-'j 

ZA,>(l-p)^pM ^ p^-p^ 

s=0 \ is=0 



which is positive as soon as j < k — 2. The quantity A( is invariant when t is replaced by 1/t, so 
we also get the desired result ior j = k — 2. 

Assume ([T4| is true for any endpoint of intervals of the subdivision ^(j), j < i — 1. Let t 
be an endpoint of an interval of -^{1); then there exists an interval [ti,t2] of J^{£ — 2) such that 
t G [ti, ^2]- We can write 

i/fe,p ([i, 00)) = Vk, p {[h, co)) + Vk^p{[ti,t2]) i^k,p{[u, 00)) 
and z/fe,p ([0, 1/t]) = i^k^p ([0, IA2]) + >^k^p ([l/t2, 1/ti]) >^k^p ([0, 1/u]) 

for some endpoint u of an interval of ^{2). If i/^^p ([ti, ^2]) > J^fe.p ([1/^2, l/^i]), we get the result 
since (|14p holds for u, and ^2- Otherwise, we can write At as 



At, - Vk^p ([ti, ts]) + Vk^p ([l/t2, lAl]) + Vk^p ([ti, t2]) i/fe,p ([m, C»)) - I/fc,p ([1/^2, 1/tl]) Vk^p ([0, 1/m]) 

which is greater than 

At, +:/fe,p([ti,t2])A„>0. 

n 
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Remark 5.2. We can also dcEne the probability measure v^p for p = I. (When k = 3, this 
is related to Minkowski's Question Mark Function, see [T].J It is straightforward to check that 
Vk,! ([i, oo)) - i^k,! ([0, 1/i]) = for all t > 0, which yields 

logxdvk.iix) = 0. 



6. Reduction: The non-linear case 

In the non-linear case, where Fn+2 ~ l-^^n+i i Fn\, the sequence {Fn)n>i can also be coded 
by the sequence (Xn)n>3 of i.i.d. random variables taking values in the alphabet {R, L} with 
probability {p, 1 — p) . Each R corresponds to choosing the + sign and can be interpreted as the 
right multiplication of (Fn-i, Fn) by the matrix R defined in ^. Each L corresponds to choosing 
the — sign but the interpretation in terms of matrices is slighty different, since we have to take into 
account the absolute value: Xn+i ~ L corresponds either to the right multiplication of (i^„_i, Fn) 
by L if (Fn-i, Fn)L has nonnegative entries, or to the multiplication by 

(., ,:-{l \ 

Observe that for all < j < fc — 2, the matrix RU has nonnegative entries (see ([5])), whereas 

RL^^^ ^ ( n 1 ) ■ Therefore, if Xi ~ R \s followed by some i's, we interpret the first (fc — 2) 

L's as the right multiplication by the matrix L, whereas the (fc — l)-th L corresponds to the 
multiplication by L' . Moreover, RL^^^L' = Id, so we can remove all patterns RL^^^ in the 
process (X„). 

We thus associate to X3 . . . x„ the word Red(a;3 . . . a;„), which is obtained by the same reduction 
as Red(x3 . . .a;„), except that the letter added in Step 1 is always Xi. We have 

(^„_i,F„) = (^i,^2)Ited(x3...x„). 

Since the reduction process is even easier in the non-linear case, we will not give all the details 
but only insist on the differences with the linear case. The first difference is that the survival 
probability of an R is positive only iip > 1/fc. 

Lemma 6.1. For p > 1/fc, the number of R's in Red(X3 . . . Xn) satisfies 

\Rcd{X^ . . . Xn)\R > +00 a.s. 

n — >OD 

and the survival probability pn is for p < 1 the unique solution in ]0, 1] of 
(16) g{x) = 0, where g{x) := (1 - x) fl + --^— x\ - 1 . 

Ifp < 1/fc, PR = 0. 

Proof. Since each deletion of an R goes with the deletion of (fc — 1) L's, ii p > 1/fc, the law of 
large numbers ensures that the number of remaining i?'s goes to infinity, li p < 1/fc, there only 
remains L's, so pn ~ 0. 

Doing the same computations as in Section 12.11 we obtain that, for all < j < fc — 2, the 
probability pj for an R to be followed by L^ R after the subsequent steps of the reduction is 

^ (1 -pYppr 

^' {l-p + PPBy+^- 

Since pn ~ X]j=o Pj' ^*^ S*^* Ih&l pji is solution of g{x) = 0. Observe that g(0) = 0, g{l) = —1, 
^'(0) > for p > 1/fc and g' vanishes at most once on R+. Hence, for p > 1/fc, pn is the unique 
solution of g{x) = in ]0, 1]. For p = l/fc, g'{0) — and the unique nonnegative solution is 
PR = 0. □ 
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6.1. Case p > \/k. As in the linear case, the sequence of surviving letters 

(5,),>3= lim Ited(X3...X„) 

n — 'oo 

is well defined for p > 1/fc, and can be written as the concatenation of a certain number s > of 
starting L's and of blocks: 

S1S2 ■ ■ ■ = L B1B2 ■ ■ ■ 
where for all i > I, Be G {R, RL, . . . , RL^^^}. These blocks appear with the same distribution 
Pp as in the linear case, but with a different parameter p. 

Lemma 6.2. In the non-linear case, for p > 1/k, the blocks {Be)i>i are i.i.d. with common 
distribution law Pp dcGned by ([9|), where p :~ '°^\/l — pr and p^ is given by Lemma \6.1\ 

As in Section 14.11 we can embed the sequence {Xn)n>3 in a doubly-infinite i.i.d. sequence 
(X*)„ga with Xn = A"* for all n > 3. We define the reduction of (Ar*)_oo<j<n by considering the 
successive Red(Ar„_Ar . . . A'„). The analog of Proposition 14. II is easier to prove than in the linear 
case since the deletion of a pattern RL^^^ does not affect the next letter. The end of the proof is 
similar. 

6.2. Case p < 1/fc. Since in this case the survival probability of an R is pr = 0, the reduced 
sequence Red(A'^) contains only L's. We consider the subsequence (Fn^) where Uj is the time 
when the j-th L is appended to the reduced sequence. This subsequence satisfies, for any j, 
-Frij+i = l^^nj — Fnj_i\, which corresponds to the non-linear case for p = 0. 

Therefore, we first concentrate on the deterministic sequence i^„+i = |Ai^„ — i^„_i|, with given 
nonnegative initial values Fq and Fi . 

Proposition 6.3. For any choice of Fq > and Fi > 0, the sequence defined inductively by 

Fn+i = \XFn - F„_i| is bounded. 

Lemma lHTSl in the next section gives a proof of this proposition for the specific case A = 2 cos Tv/k. 
We give here another proof based on a geometrical interpretation, which can be applied for any 
0< A<2. 

The key argument relies on the following observation: Let be such that A = 2 cos 0. Fix two 
points Po, Pi on a circle centered at the origin O, such that the oriented angle (OPq, OPi) equals 
6. Let P2 be the image of Pi by the rotation of angle 9 and center O. Then the respective abscissae 
xo, xi and X2 of Pq, Pi and P2 satisfy X2 = Xxi — xq. We can then geometrically interpret the 
sequence (P„) as the successive abscissae of points in the plane. 

Lemma 6.4 (Existence of the circle). Let e]0, 7r[. For any choice of [x, x') G R^ \ {(0, 0)}, their 
exist a unique R > and two points M and M' , with respective abscissae x and x' , lying on the 
circle with radius R centered at the origin, such that the oriented angle {OM, OM') equals 9. 

Proof. Assume that x > 0. We have to show the existence of a unique R and a unique t € 
] — 7r/2,7r/2[ (which represents the argument of M) such that 

Rcost = X and R cos{t + 9) — x' . 

This is equivalent to 

x' 
R cos t = x and cos 9 — tan t sin 9 = — , 

x 

which obviously has a unique solution since sin 9 ^ 0. 

If a; = 0, the unique solution is clearly R = x' / cos{9 — tt/2) and t = —tt/2. 

Remark: Since xi > 0, we have t + 9 < 7r/2. D 



Proof of Provosition 1 6. 3i At step n, we interpret P„+i in the following way: Applying the lemma 
with X = P„_i and x' = P„, we find a circle of radius P„ > centered at the origin and two 
points M and M' on this circle with abscissae x and x' . Consider the image of A/' by the rotation 
of angle 9 and center O. If its abscissa is nonnegative, it is equal to Pn+i, and we will have 
Rn+i = Rn- Otherwise, we have to apply also the symmetry with respect to the origin to get a 
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Ti/k 




Fn+2 Fn+l Fn~l Fn 

Figure 2. Rn — Rn+i is the radius of the largest circle, and Rn+2 is the radius 
of the smallest. 



point with abscissa Fn+i- The circle at step n + 1 may then have a different radius, but we now 
show that the radius always decreases (see Figure [2|). 

Indeed, denoting by a the argument of M', we have in the latter case tt/2 — 9 < a < tt/2, 
Fn = Rn cos a and Fn+i = i?„ cos(a + 9 + n) > 0. At step n + 1, we apply the lemma with 
X = Rn cos a and x' = Rn cos(a + + 7r). From the proof of the lemma, if Fn = (i.e. if a = 7r/2), 
Rn+i = i?„cos(a + 6' + tt)/ cos{9 — n/2) = Rn- If Fn > 0, we have i?„+i ~ Rncosa/ cost, where 
t is given by 

cos(a + 9 + tt) 



cost 



tant siii9 = 



(cos 9 — tan a sin ( 



cos a 
We deduce from the preceding formula that tant + tana — 2 cos 9/ sin 9 > 0, which implies t > —a. 
On the other hand, as noticed at the end of the proof of the preceding lemma, t + 9 < 7r/2, hence 
t < a. Therefore, cos a < cost and Rn+i < Rn- 

Since Fn < Rn < Ri for all n, the proposition is proved. D 



IXFn — Fn-i\ and its two first 



We come back to the specific case A = 2cos7r/fc. 

Proposition 6.5. Let (Fn) be inductively defined by Fn+i 
positive terms. Tfie following properties are equivalent: 

(1) Fq/Fi admits a finite X-continued fraction expansion. 

(2) The sequence {Fn) is ultimately periodic - 

(3) There exists n such that Fn = 0- 

Proof- We easily see from the proof of Proposition 16.31 that (2) and (3) are equivalent. We now 
prove that (3) implies (1) by induction on the smallest n such that F„ = 0. If F2 ~ 0, then 
|AFi — Fo\ = 0, and we get Fq/Fi = A. Let n > 2 be the smallest n such that F„ = 0. By the 
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induction hypothesis, F1/F2 admits a finite A-continued fraction expansion. Therefore, 



Fi F1/F2 

admits a finite A-continued fraction expansion. 

It remains to prove that (1) impUes (3). We know from Proposition 13.31 that aU positive real 
numbers that admit a finite A-continued fraction expansion are endpoints of generalized Stern- 
Brocot intervals, hence by (jlOp . can be written as [1, ai, . . . , ojJa with a^ = ±1 for any i and 
such that we never see more than [k ~ 1) alternated ±1 in a row. We call such an expansion a 
standard expansion. Conversely, all real numbers that admit a standard expansion are endpoints 
of generalized Stern-Brocot intervals, hence are nonnegative. Assume (1) is true. If i^o/^i = [1]a, 
then F2 = 0. Otherwise, let [1, ai, . . . , aj\\ be a standard expansion of Fq/Fi. Then, 

^1 1 r 1 

F2 \X-Fo/F,\ ^ ' ' ^J 

If ai = 1, then [oi, . . . , aj]x > and it is equal to F1/F2. Otherwise, F1/F2 = [— ai, —02, . . . , —aj]x. 
In both cases, we obtain a standard expansion of F1/F2 of smaller size. The result is proved by 
induction on j . D 

Remark 6.6. In general, if Fq/Fi does not admit a Unite X-continucd fraction expansion, (F„) 
decreases exponentially fast to 0. However, the exponent depends on the ratio Fq/Fi. 

We exhibit two examples of such behavior. 

Let q := {X + VA2 + 4)/2 be the fixed point of /o. Start with h/Fi = q. Then, by a 
straightforward induction, we get that for all n > 0, F„ = q^"-Fo. 

Start now with Fq/Fi = q' , where q' is the fixed point of /i. Then, we easily get that for all 
" > 0, F2n = [q' h{q'))~'^FQ and F2n+i = F2n/q'- The exp onent is thus I / ^JYMff) , which is 
different from \/q: For k ^ 3, q ^ (j) (the golden ratio) and \/q' fo{q') = \/0- 

Proof of Theorem \1.2[ We have seen that the subsequence (Fn), where Uj is the time when the 
j'-th L is appended to the reduced sequence, satisfies, Fn-^-^ = \XFn- — Fn_-^\ for any j. From 
Proposition l6.3[ this subsequence is bounded. Moreover, we can write Uj ~ j + kdj, where dj is the 
number of _R's up to time Uj. By the law of large numbers, dj/rij -^ p, and we get j /nj ^ 1 — kp. 
This achieves the proof of Theorem 11.21 D 

7. Case A > 2 

The case A > 2 (p > 0) is even easier to study since there is no reduction process. 

Observe that the linear and the non-linear case are essentially the same. Indeed, in the non- 
linear case, F(F„+i/F„ > l|F„_i,^„) > p and if F„+i/F„ > 1, then Fn+2/Fn+i > T Therefore, 
with probability 1, there exists 7V+ such that for all n > N^, the quotients F„+i/F„ are larger 
than 1. Moreover, for n > N^, there is no need to take the absolute value and the sequence 
behaves like in the linear case. We thus concentrate on the linear case. 

We now fix A > 2. The sequence of quotients Qn '■— Fn/ Fn-i is a real-valued Markov chain 
with probability transitions 



Qn = q] = 1 -p, 



P [Qn+i = fniq) Qn=q) =P and P f Q„+i = fL{q) 
where fniq) := X + l/q a nd /^(q) := A - 1/q. 

Let B := G [1, A] be the largest fixed point of /i. Note that we have P((5„+i > 

X\Qn) > min(p, 1— p) for any n > 2 and, again, if Qn > B, then Qn+i > B. Thus, with probability 
1, there exists N+ such that for all n > N+, the quotients Qn are larger than B. Without loss of 
generality, we can henceforth assume that the initial values a and b are such that Q2 > B. 
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Wc inductively define sub-intervals of R-|_ indexed by finite sequences of i?'s and Us: 

1 



'R 



fR{[B,^]) 



A,A + 



B 



and /l :=/l([S,oo]) = [S,A], 



and for any finite sequence X in {R,L}*, 

IxR ■= Ir{Ix) and Ixl := Il{Ix)- 
Obviously, all these intervals are included in \B, A + -g] . 

Lemma 7.1. Let W and W he two Unite words in {R,L}*. 

• IfW is a suffix ofW, then Iw C Iw; 

• If neither W is a suffix of W nor W is a suffix of W, then Iw and Iw' have disjoint 
interiors. 

Proof. The first assertion is an easy consequence of the definition of Iw- To prove the second 
one, consider the largest common sufHx S oi W and W . Since LS and RS are suffix of W and 
W, by the first assertion, it is enough to prove that Ils and Irs have disjoint interiors. This 
can be shown by induction on the length of S, using the fact that fa and /^ are monotonic on 

[B,oo]. D 

Lemma 7.2. Let {Wi)i>i be a sequence ofR's and L's. Then C\n>i Iw„...Wi j's reduced to a single 
point. 



Proof. By Lemma 17.11 /w„+iW„...Wi C Iw^.-Wi- Since the intervals are compact and nonempty, 
their intersection is nonempty. It remains to prove that their length goes to zero. First consider 
the case A > 2. The derivatives of /l and //? are of modulus less than 1/B^ < 1. Therefore, 
the length of Iw„...Wi is less than a constant times (1/i?^)". Let us turn to the case A = 2. 



j+i 



' 3 



which is of length 1/j. Hence, if W, 
^Wi is included, for some r < n, in luWr-.-Wx = Iwi ° 



Wi contains j consecutive 
o fw^ [Ihi ) which is of 



Observe that /^. 

L's, then Iwn...Wi ^^ mciuueu, iui Bomt; / ^ /t, in J^L^Wr-.-W 

length less than l/j (recall that the derivatives of /l and /^ are of modulus less than 1). On the 
other hand, the derivatives of /l o fp and /^ o ffj are of modulus less than 1/{2B + 1)^ = 1/9 on 
[B, oo]. Therefore, considering the maximum number of consecutive L's in Wn ■ ■ ■ Wi, we obtain 

SUPvy„...Wi \Iw„...Wi\ >■ 0. D 



4/3 



Illl ' Ir 



3/2 

7/5 11/7 



{1-pf (l-pfp 



III 



(l-p)' 



8/5 



5/3 



II 



Irrl Ilrl 



(i-p)p' a-pfp 



Irl 



{l-p)p 



1 — p 



7/3 



12/5 



5/2 

17/7 13/5 



8/3 



^LRR Ir 



Irr 



Irl r Illr 



Ilr 



{l-p)p 



Ir 



Figure 3. First stages of the construction of the measure /J.p_2- 
We deduce from the preceding results the invariant measure of the Markov chain (Qn). 



18 ELISE JANVRESSE, BENOIT RITTAUD, THIERRY DE LA RUE 

Corollary 7.3. The unique invariant probability measure fipx of the Markov chain (Qn) — 
(Fn/Fn-i) is given by 

(17) Mp,A(/w):=pl^l''(l-p)l'^l^ 

for any finite word W in {R,L}*, where \W\ii and \W\l respectively denote the number of R's 
and L's in W. 



We can now conclude the proof of Theorem 11.31 by invoking a classical theorem about law of 
large numbers for Markov chain (see e.g. [B], Theorem 17.0.1). 

Note that the explicit form of the invariant measure when p = 1/2 and A > 2 was already given 
by Sire and Krapivsky |10| . 

8. Variations of the Lyapunov exponents 

8.1. Variations with p. 

Theorem 8.1. For any integer fc > 3, the function p i— > ^p.x^ is increasing and analytic on ]l/fc, 1[, 
and the function p ^-^■ ^p.x^. is increasing and analytic on ]0, 1[. Moreover, 

(18) lini7p,A, = lim >,A, =0, 

and 



(19) lim 7p,A,. = liM = 1™ 7p,Afc = 7i,a,. = log — ^ . 

P^l p^l \ Z J 

For any A > 2, the function p i— > 7p^A is increasing and analytic on ]0, 1[. 

The proof of the theorem relies on the following proposition, whose proof is postponed to the 
end of the section. 

Proposition 8.2. Let (Xi) be a sequence of letters in the alphabet {R, L} and {XD be a sequence 
of letters in the alphabet {R,L} obtained from (Xi) by turning an L into an R. If X = A^ for 
some k > 3, then, in the non-linear case, any label Fn coded by the sequence (Xi) is smaller than 
the corresponding label F^ coded by {X[). If X>2, and if F2/F1 > 1, any label Fn coded by the 
sequence (Xi) is smaller than the corresponding label F^ coded by {X'^). 



Proof of Theorem \8.1[ Let A = Afe for some integer fc > 3. Let 1/k < p < p' < 1. Let (Xi) 
(respectively {X^)) be a sequence of i.i.d. random variables taking values in the alphabet {R,L} 
with probability {p, I — p) (respectively {p' ,1 — p')). We can reahze a couphng of {Xi) and {X'^) 
such that for any i, Xi = R implies X'i ~ R. From Proposition 18.21 it follows that the label F„ 
coded by {Xi) is always smaller than the label F^ coded by {X'^). We get that 

7p,Afc = hm - log Fn < hm - log F^ = %,^x, ■ 

Therefore, p 1— > 7p,Afc is a non-decreasing function on [1/fc, 1]. 

Observe that p h^ pji is non-decreasing in both (linear and non-linear) cases. Hence, the 
function p : p t-^ ''^\/l — Pr is non-increasing in both cases. We conclude that p i-^ lp,\k is 
non-decreasing on [0, 1]. 

Since 7p,Afc > for < p < 1, the upper Lyapunov exponent associated to the product of 
random matrices is simple, and we know from [7j that 7p.Afc is an analytic function of p g]0, 1[, 
thus it is increasing. Via the dependence on p which is an analytic function of p, we get that 7p.Afc 
is an analytic increasing function of p e]l/fc, 1[. 

Now, observe that p \ — > J„ \ogxdvk.p{x) is continuous on [0, 1] (as the uniform limit of con- 
tinuous functions). When p goes to zero in the linear case (or p -^ 1/fc in the non-linear case), 
pn tends to and p tends to 1. By continuity of the integral, we obtain (|18|1 using Remark 15.21 
When p = 1, the deterministic sequence Fn ~ Fn grows exponentially fast, and the expression of 
7i^Afc follows from elementary analysis. 
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When A > 2 (wc do not need to distinguish the hncar case from the non-linear cases) , the proof 
is handled in the same way, using Proposition 18.21 D 



Proof of Proposition \8.8\ when A > 2. We let the reader check that in this case, for all s > the 
matrix RL'^ has nonnegative entries. Suppose the difference between (Xi) and (Xl) occurs at 
level j. For any n > j, the sequence Xj . . . Xn can be decomposed into blocks of the form RL^, 
s > 0, hence the product of matrices Xj ■ ■ ■ X„ has nonnegative entries. If n > j, we can thus 
write F'^ as a linear combination with nonnegative coefficients: F^ = CiF',_2 +C2F',^. Moreover, 
F„ = -CiFj_2 + Cai^j-i = -Ci^i-2 + C2F,'_i, hence F„ < F^^ (since F2/F1 > 1, all F„'s are 
positive) . n 



The proof of Proposition 18.21 when A = Afc uses three lemmas. The first one can be viewed as 
a particular case when the sequence of i?'s and L's is reduced. 

Lemma 8.3. Let A = Afc. Let a > 0, b > 0, ji > and J2 > such that ji + I + J2 < k — 2. If 
{a',b') ^ {a,b)RL3^RL3^ and {a" ,b") = {a,b)RU^+^+^\ then b' > b" . 

Proof For any €e {0,...,j2}, set (x,,a;f+i) := {a,b)RU^RL\ and [yi^yt+i) := {a,b)RU^+^+'. 
Then the quotient x^+i/xf lies in li (see Section|3]), whereas the quotient y^+i/y^ lies in /jj+i+£. It 
follows that yi^i/yi < xi+i/xi^ and since xq = yo, we inductively get that for all £ G {0, . . . , J2 + I}, 
yi < Xi. The lemma is proved, observing that b' = Xj^+i and 6" = yj-^+i- D 

Lemma 8.4. Let A = Afc. Let {Xi)i>2 be a sequence of matrices in {R,L}, which does not 
contain fc — 1 consecutive L's and such that X2 = R. Let xq > 0, xi > 0, and set inductively 
{xi,x,+i) := (xj_i,Xi)X,+i. Then for any i > 0, Xi+k > Xi. 

Proof. If Xi^i = R, this is just a repeated application of the following claim: If a > 0, 6 > 0, 
< j < fc - 3, and if we set (a', 6') := {a,b)RL^ , then b' > b. Indeed, by ©, wc have b' > 
b sm{{j + 2)7r/fc)/ sin(7r/fc) > b. 

If Xi^i ~ L, wc first prove the lemma when the sequence Xi^i . . . X^+fc contains only one 
R: Xi+j = R for some j e {2,...,fc — 1}. We proceed by induction on j. If j = 2, then 

(xi+fc_i, Xi+fc) — (xi^i, Xi)LRL''^'^ . By ([5]), the second column oi RL^~^ is ( „ ), thus x^+fc — Xi. 

Now, assume j > 2 and that we have proved the inequality up to j — 1. Since the sequence 
of matrices starts with an R and docs not contain fc — 1 consecutive i's, we have x^+i/xi G li 
for some i < k ~ j (see Section [3]). In particular, x^+i/x^ > b^-j- Now define x'^.j^^^ by 
ixi+k,x[^k+i) '■= ixi+k-i,Xi+k)L. We have x[^^._^_i/xi+k e h-j+i, thus is bounded below by 
bk-j- Using the induction hypothesis x'^^j,^^ > Xi+i, we conclude that x^+fc > x,;. 

Finally, assume that the sequence X^+i . . . Xi^k starts with an L and contains several i?'s. 
Turning the last R into an L, we can apply Lemma. 18.31 to compare x^-i-fc with the case where there 
is one less R, and prove the result by induction on the number of i?'s. D 

Lemma 8.5. Let A = Afc. Let Fn be inductively defined by Fq > 0, Fi > and Fn+i ~ 
\XFn — Fn^\\ for any n > 1. Then for any n > 0, Fnj^k < Fn- 

Proof. For n < 0, let G„ := F_„ > 0. Then, for any n < —1, we have 

(T* r* \ — J (^"-1'^")^ if AF„ > F„_i, 
]{Gn-i,Gn)R otherwise. 

Moreover, we can assume that the sequence of matrices in {R, L} corresponding to (G„) never 

contains fc — 1 consecutive L's. Indeed, the second column of i'"'^^ is I „ ]. Thus, if we had 

fc — 1 consecutive L's, we could find n such that —Gn-i — Gn+k-i, which is possible only if 
Gn-i = Gn+k-i = 0. But if such a situation occurs we can always turn the first L into an 
R without changing the sequence (because (0, G„)i? = (0,G„)L). The result is thus a direct 
application of Lemma 18.41 D 
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Proof of Provosition \8.'t^ when A = Afc. Suppose the difFerence between (Xi) and (X,') occurs at 
level j. We decompose {Xj)i>j as LL^Y and {X'Ai>j as RL'^Y, where < r < +00 and Y — 
{Yi)i>j+r+i is a sequence of letters in the alphabet {R,L} such that Ij+^+i = R- 

Suppose first that, after the difference, all letters are L's {Y = 0). Let ji € {0, . . . , fc — 2} be 
such that Fj_i/Fj e Ij^. Without loss of generality, wc can assume that the sequences {Xi) and 
{X'i) are reduced before their first difference. Then, Xj-j^-i . . . Xj^i = X',_,_^ . . . X';-^ = RU^ . 

By Lemma [5751 Fjj^s > ^j+s fo'^ ^^1 < s < j2, where 72 := A: — 3 — ji. 

Now, by Lemma [8.41 for all 1 + J2 < s < fc — 2, Fj+s > Fj+s-k, which is equal to F'-_^^^_f, 
since s < k. On the other hand, when 5=^2 + 1, we have F'-^-^^_j^ = ^,'+,,+1 because 
X^_j^_i . . . Xj+j,+i = RL''-\ Moreover, by LemmaJHH F^+,_k > ^j+, for alll+ja < s < fc-2. 
We thus get that Fj+s > F^+s for all J2 + 1 < s < fc - 2. 

If s > fc — 1, reducing the pattern RL'^^^ in the sequence {Xj)i>j, we have -Fj+s = Pj+s-k 
which is larger than F'^^ by Lemma 18.51 

Suppose now that the suffix Y is reduced. The above argument shows that all labels up to 
j + r are well-ordered: In particular, Fj_^.j.-i < F'-_^^^_^ and -Fj+r ^ ^j+r- Since Y is reduced, 
we can write, for any n > j + r, (_?„, -Fn+i) = (-Fj+r-i, Fj+r)Yj^r+i • ■ ■ ^n+i, where each Yi is 
interpreted as the corresponding matrix (the same equality is valid if we replace F by F'). The 
product Yj^r+i ■ ■ -Yn+i can be decomposed into blocks of the form RL^, with < ^ < fc — 2, 
which are matrices with nonnegative entries. Therefore, for any n > j + r, the label Fn is a linear 
combination of i^j+^-i and -Fj+r , with nonnegative coefficients. Moreover, it is also true with the 
same coefficients if we replace F hy F' . We conclude that Fn < F!^. 

In the general case, we make all possible reductions on Y . Wc are left either with a reduced 
sequence or with a sequence of L's, which are the two situations wc have already studied. D 

Remark 8.6. In [5], a formula for tie derivative ofjp^i with respect to p was given, involving the 
product measure v^^p ^ 1^3, p- We do not know whether this formula can be generalized to other 
fc's. 

8.2. Variations with A. Forp = 1, the deterministic sequence F„ ~ Fn grows exponentially fast, 
and wc have in that case 



. /A + ^/A^T4^ 
7l,A == 7i,A = fog I ^ I , 

which is increasing with A. 

We conjecture that, when p is fixed, jp^x^. and 7p,Afc are increasing with fc, and that 7^^^ is 
increasing with A for A > 2 (see Figure |3|). 

9. Connections with Embree-Trefethen's paper 

9.1. Positivity of the Lyapunov exponent. We have proved that the largest Lyapunov expo- 
nent corresponding to the linear A-random Fibonacci sequence is positive for all p. In [5] , Embree 
and Trefethen study a slight modification of our linear random Fibonacci sequence when p = 1/2. 
To be exact, they study the random sequence Xn+i = a;„ ± (3xn-i, which by a simple rescaling 
gives our linear A-random Fibonacci sequence where A = 1/%//? (see our introduction). However, 
the exponential growth is not preserved by this rescaling. More precisely, the exponential growth 
cr(/3) = lim |a;„|^/" of Embree and Trefethen's sequence satisfies 

logcr(/3) =7i/2,A-fogA. 

In particular, cr(/3) < 1 if and only if 71/2, a < fog'^i which according to the simulations described 
in their paper happens for /? < /3* ss 0.70258 . . . (which corresponds to A > 1.19303 . . .). 

By Theorem 18.11 the function p 1— > 7p,\ is continuous and increasing from to 71. a > log A. 
Hence there exists a unique p*{X) G [0, 1] such that, for p < p* , ^p^\ < log A and for p > p* , 
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7p,A 




Figure 4. The value of ^p^x (linear case, left) and 7p.A (non-linear case, right) 
for A = Afc, fc = 3,4,5,10, A = 2 (bold), A = 2.05, 2.1, 2.5 and 3. Numerical 
computations support the conjecture that ^p^\ and ^p,\ are increasing with A. 



7p,A > log A. According to [2], for A = 1 we have p* < 1/2, and for A = Afc {k > 4) and A > 2, 
p*'> 1/2. 

For A > 2, we can indeed prove that 71/2, a < log A: By Jensen's inequality, we have 

/ .A+l/B \ 

7i/2,A < log I / xdfii/2,\ 1 , 

which is equal to log A by symmetry of the measure ^1/2. a ■ 

For A = 1, we know that ■jp^i > for all p > thus p* = 0. When A = Afc, k > 4, numerical 
computations of the integral confirm that p* > 1/2, but we do not know how to prove it. 

9.2. Sign-flip frequency. Embree and Trefethen introduce the sign-flip frequency as the propor- 
tion of values n such that FnF„+i < 0, and give (without proof) the estimate 2^'^'^/^**^'^ for this 
frequency, as A ^ 2, A < 2. 

Note that, for A > 2, there are no sign change as soon as n is large enough, and the sign-flip 
frequency is zero. 

For A = Afc, recall that for n large enough, the sign of the reduced sequence {F^) is constant 
(see Lemma [3. 2p . Moreover, by ^ and the fact that for all < j < A: — 2 the matrix RV has 
nonnegative entries (see ^), the product FnF^ changes sign if and only if a pattern RL^^^ is 
removed. Thus, the sign-flip frequency is equal to the frequency of deletions in the reduction 
process. 

Note that we have to make sure that this frequency indeed exists. This can be seen by consider- 
ing the reduction of the left-infinite i.i.d. sequence {X*)^_^ (Section l4.ip . since for n large enough, 
deletions in the reduction process of {X)"^ occur at the same times as in the reduction process of 
(X*)"q^. In the latter case, the ergodic theorem ensures that the frequency a of deletions exists 
and is equal to the probabihty that [X*)^_^ be not proper. By Lemma H?^ {X*)'^_^ is not proper 
if and only if there exists a unique ^ > such that {X*)'^_^ is an excursion, and {X*)Zo^^ is proper. 
Thus, 

^^ P(u;)(l-cr). 



E 



w excursions 



By ()lip . we get that the sign- flip frequency is equal to 

P(l -Pr) 



(20) 



c^(Afc,p) 



p + (1 - p)pB. + p{l - PR.) ' 

Now, for a fixed p g]0, 1[, we would like to obtain an estimate for cr as fc ^- cx3. First, observe 
that PR = PR(k) -^ 1 as k ^> <x). Indeed, recalling the expression of the function g given by ([S]), 
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for any x g]0, 1[, wc have g{x) < for k large enough, which imphes pn > x. Then, since pn 

satisfies 

/ \ fc-i 

1- = ( 1- PPfl 

we get that p^i — > 1 exponentially fast with k. Using this estimation in the above equation, 
elementary computations lead to 

Thus, 

a(Afe,p)^ p{l-pf-\ 

k — >oo 

For p = 1/2, this proves the estimate provided in [2] in the special case A = Afc. 
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